Finite temperature properties of the two-dimensional SU(2) Kondo-necklace 
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We analyse several thermodynamic properties of the two-dimensional Kondo necklace using finite- 
temperature stochastic series expansion. In agreement with previous zero-temperature findings the 
model is shown to exhibit a quantum critical point (QCP), separating an antiferromagnetic from a 
paramagnetic dimerized state at a critical Kondo exchange-coupling strength J c ~ 1.4. We evaluate 
the temperature dependent uniform and staggered structure factors as well as the uniform and 
staggered susceptibilities and the local 'impurity' susceptibility close to the QCP as well as in the 
ordered and quantum disordered phase. The crossover between the classical, renormalized classical, 
and quantum critical regime is analyzed as a function of temperature and Kondo coupling. 
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There is growing evidence that unconventional finite 
temperature properties of many novel materials stem 
from zero-temperature phase transitions, i.e. changes of 
the ground state symmetry as a function of some control 
parameter. Prominent potential candidates to show such 
quantum phase transitions are the cuprate superconduc- 
tors Q, Q , quantum magnets 0, fj| , and heavy-fermion 
or dense Kondo systems 0, 0]. In the latter, quantum 
critical points (QCPs) can arise from the competition be- 
tween magnetic long-rage ordered (LRO) and renormal- 
ized paramagnetic metallic or semimetallic phases result- 
ing from local Kondo-screening. This has been conjec- 
tured early on by Doniach Q. Semimetallic behavior in 
nonmagnetic states of Kondo lattice materials is typical 
for stoichiometric "Kondo insulators" like CeNii-^Pt^Sn 
which undergoes a para-to-antiferromagnetic transition 
at x w 0.2. ..0.3 [3, LJ|- A model for such materials is 
Kondo-Hubbard lattice model (KHLM) 



■jJ^SpiSn (!) 

l,a/3 



with conduction electrons c la of spin S pi , which are cor- 
related via and on-site Coulomb repulsion U, and coupled 
by antiferromagnetic (AFM) Kondo-exchange to local- 
ized spins Sn at sites I. At half filling on bipartite lat- 
tices in D> 2 dimensions and in the strong-coupling limit 
U/t > 1 the HKLM shows AFM LRO if the conduction- 
electron superexchange j ~ t 2 /U dominates the Kondo 
scale Tk ~ t exp (— 1/pJ) where p refers to the DOS 
[Icj . Kondo screening will prevail if Tx/j ^S> 1. On 2D 
square lattices the critical coupling j c (U/t) has been de- 
termined at temperature T = using projector QMC 
and bond-operator methods 0|. In the strong cou- 
pling limit and at half filling Eqn. (TIJ simplifies to the 
SU (2)-symmetric so-called Kondo-necklace (SKN) 

Hskn = j ^ Sp; • Sp m + jJ2 s Pi-Sn (2) 

lm I 

with j = 1 hereafter. In this work we will focus on the 2D 
square lattice, where at T = the SKN has been inves- 
tigated by bond-operator methods, series expansion and 
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Figure 1: Solid: staggered structure factor S n (Q) vs. J close 
to the QCP at low-temperatures and for two system sizes: a) 
T = 0.05, L = 24, b) T = 0.01, L = 24, and c) T = 0.05, 
L = 34. Dashed: fits of S P (Q) to c [J c - J] v for J > 1 with 
J c indicated per panel. The fits depart visibly from Sp (Q) 
only close to J w 1.4. Dashed-dotted: ratio of total staggered 
structure factor to susceptibility times T. If not indicated, 
statistical errors are less than the solid-circle marker size. 



exact diagonalization H3. Il5|. These studies located 
a QCP at J c - 1.370 . . . 1.408 separating AFM LRO from 
a gapped spin-dimer phase. The latter can be viewed as 
the strong-coupling analog of the Kondo-screened para- 
magnetic state of the KHLM. 

While the ground state properties of the 2D SKN have 
been studied rather extensively, thermodynamic and fi- 
nite temperature critical properties of the SKN remain an 
open issue. Therefore the aim of this work is to shed light 
on the 2D SKN at finite temperatures using a quantum 
Monte-Carlo (QMC) approach. To this end we employ 
the stochastic series expansion (SSE) with loo p-upda tes 
introduced by Sandvik and Syljuasen in Refs. [la, 03 to 
which we refer the reader for details on this approach. 

We start by discussing the longitudinal staggered 
structure factor 



Sn (Q) 



(3) 



where 



= J]; exp (iQ • r/) /iV„ is the staggered 
magnetization with Q = (ir,ir,ir). m^q selects between 
n = P,I,A, for which 17 runs over the 'conduction elec- 
tron plane' for n = P, the 'Kondo sites' sites for n = I, 
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Figure 2: Uniform susceptibility Xu vs - temperature for 
0.05 < T < 10 and 0.5 < J < 2.0 at L =24 (solid) and 
for 1.39 < J < 1.42 at L = 50 (100) (dashed(dashed-dotted)). 
Statistical errors are less than the solid-circle marker size. 
The difference between L =50 and 100 at J = 1.40 remains 
below the statistical error for all T depicted. Legends label 
plots from top to bottom 

and all sites for n = A. Fig. ^ shows the squared stag- 
gered moment Mq = 3S n (Q) vs. J at low temperatures. 
The system sizes Na are L x L x 2 = N with periodic 
boundary conditions (PBC) in the planar directions and 
L = 24 and 34. In all three panels Mq is finite below a 
critical value of J — J c and drops to approximately zero 
for J > J c . We identify J c with the QCP and expect 
AFM LRO for J < J c in the thermodynamic limit at 
T = 0. For J > J c we find no other transitions, i.e. the 
systems connects adiabatically to the limit of J = oo. 
Therefore it is dimerized. At fixed N, Mq will saturate 
for T — > due to finite size gaps. For L = 24 this is the 
case in Fig. \l\b) for J > 1. Fig. fallows no conclu- 
sion about the magnitude of the T = order parameter, 
which requires finite size scaling analysis yjj- The crit- 
ical coupling however can be extracted efficiently from 
these results since J c is almost invariant to increasing N 
or lowering T relative to the parameters in Fig. To 
determine J c we fit Mq to a power law Mq ss c [ J c — J]' y 
in its region of negative curvature and for J > 1. This 
procedure depends only little on the interval of J fitted 
to. The resulting scatter of J c is taken to be a mea- 
sure of the error and is displayed in Fig. Ela)-c). We 
find that J c « 1.41(2). This agrees with J c « 1.41(1.39) 
from T — series-expansion [13] (0|) an d is also close to 
J c » 1.37 from bond-operator Bruckner theory [LH Hg|. 
A critical value of J c « 1.4 (At 2 /U) is also consistent with 
projector-QMC at T = for the KHLM [HI- 

For small J, S n (Q) remains strongly temperature de- 
pendent down to T <C 1 which is due to the near de- 
coupling of the /-sites from the planar sites leading to a 
Curie-like contribution which is cutoff only at very low T. 
This is visible already at J < 0.6, by comparing Si (Q) 
in panels a) and b) of Fig. ^ 
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Figure 3: Low-temperature uniform susceptibility X" CO 
close to the QCP for J = 1.39, 1.40, and 1.42 (from top to 
bottom). Dotted curves: fits of \u (T) to a + bT + cT 2 . Sys- 
tem sizes L =24 (solid), 50 (dashed), 100 (dashed-dotted) are 
indicated for T <0.1. Statistical errors are less than the solid 
circle marker size. For T > 0.11 finite size effects are below 
statistical error and results for J = 1.40(1.39, 1.42) refer to 
L = 100(24) only. The difference between L =50 and 100 at 
J = 1.40 remains below the statistical error for all T depicted. 

In addition to S n (Q) Fig. includes results for the 
longitudinal staggered susceptibility 



Xn(Q) 



dr (r<Q (t) 7< q ) 



(4) 



for n = A which have been encapsulated in the ratio 

Sa (Q) 



R 



T X A (Q) 



(5) 



This ratio relates the analysis of \a (Q) to that of the 



AFM non- linear a- model (NLctM) (20, |2lJ, |22|. From 
there it is expected that in the classical high-T, as well as 
in the low-T renormalized classical regime R = 1, while 
R ^ 1 in the quantum critical regime. While this is 
consistent with R (J) in Fig. we will clarify later that 
the deviations of R from unity for J « J c are strongly 
affected by finite size effects. 

Next we discuss the uniform susceptibility 



Xu = f3(m 2 ) 



(6) 



where m = J2i ^ai/^a is the total magnetization. In 
contrast to Eqn. yj, the uniform susceptibility reduces 
to a simple expectation value, since [Hskn, m] — 0. Fig. 
Elis a log- log plot of the dependence of Xu on temperature 
over more than two decades 0.05 < T < 10 and for 0.5 < 
J < 2 with system sizes L = 24, 50 and 100. For T > 0.1 
finite-size effects are negligible if L > 24. For 0.05 <T< 
0.1 finite-size effects, albeit small, have been considered 
for J in the vicinity of the QCP. As can be seen from 
the near identity of results with L — 50 and 100 at J = 
1.40 in Fig. [2J it is sufficient to choose L > 50 to reach 
the thermodynamic limit for all temperatures studied. 
For T > 2 the uniform susceptibility turns Curie-like, 
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Figure 4: Uniform susceptibility Xu vs. J 011 approaching 
the QCP for two temperatures T = 0.05 (0.01) (solid(dashed) 
line) at fixed L — 24, and comparing two system sizes L — 
24 (34) (solid(dashed-dotted) line) at fixed T = 0.05. Finite 
size effects are exceedingly small and are shown in the inset for 
better visibility. Statistical errors are less than the solid-circle 
marker size. 

independent of J. For J > J c the spin-spectrum develops 
a gap A which implies a low-temperature behavior Xu °c 
cxp (— 0A). This is consistent with Fig. [3J where to 
within statistical error Xu (T — 0.05, J = 2) = 0. For < 
J < J c AFM LRO occurs at T = 0, which agrees with 
the saturation of Xu {T — > 0) = x u shown in the figure. 
We note, that as J vanishes xl wm diverge due to the 
Curie contribution form the impurity spins. 

At the QCP we expect scaling of the uniform suscep- 
tibility. Indeed, for J « J c , and at low temperatures 
Xu follows nearly straight lines in Fig. [2J A close- 
up of this low-T region, shown in Fig. evidences a 
weak curvature of Xu{T) independent of the system size. 
These results allow excellent fits to a scale-free behav- 
ior of the form Xu ~ a + bT c with c « 1.25. However, 
this exponent differs from that obtained in the NLctM, 
i.e. c = 1 EU. Assuming the SKN to be of the 
same universality class than the NLcrM we are forced to 
treat the curvature in Fig. Elas deviations from scaling 
present already at rather low temperature. As is shown 
in Fig. 03 a reasonable description of the QMC results 
can be obtained including a second order nonuniversal 
contribution. This behavior should be contrasted against 
the AFM bilayer Heisenberg model, where critical linear 
T-scaling has been found in a comparable temperature 
range 23]. At the QCP x« vanishes for T — > due to 
the opening of the spin gap and for J > J c exponential 
behavior should replace the scaling. Vanishing of the off- 
set a at J c = 1.40(1) in Fig. Elis consistent with as from 
the static structure factor. 

For the sake of consistency it is interesting to consider 
the uniform susceptibility Xu also as a function of J in 
the low-temperature limit T <C J '. On approaching the 
QCP from the LRO side one expects Xu to vanish due to 
the incipient spin gap. The corresponding QMC results 
are shown in 31 Extracting J c from this figure is less 
straightforward, both, due to the sizeable temperature 



top to bottom: J- J 




o.i T 1 



Figure 5: Impurity susceptibility x'° c vs. temperature. Solid 
line with solid circle markers: L = 24 for 0.05 < T < 10 
and 0.5 < J < 2.0. Dashed line with solid circle markers: 
L = 50 for 0.05 < T < 0.1 and 1.39 < J < 1.42. Inset: low-T 
region including additional results for L = 40 at J = 1.40 
(Dashed-dotted). Legends refer to lines from top to bottom. 
Statistical errors are less than the solid-circle marker size in 
the main panel and are indicated by bars in the inset. In the 
quantum critical regime x'° c displays a cross-over region with 
power-law behavior x'° c ~ cT~ a with a « 0.20(5). 

variation and to the lack of a scaling prescription for Xu 
as J — > J c - Nevertheless, as can be seen in the inset, 
a value of J rj 1.4 for the QCP is consistent with the 
suppression of Xu- 

Now we turn to the individual impurity-spins longitu- 
dinal susceptibility 

X l ° c = f \r{T T Sh{T)Sh) (7) 
Jo 

where I refers to a particular site, say I = within the 
'Kondo spin' layer /. Fig. shows a log-log plot of x l °° 
vs. T for 0.05 < T < 10 and 0.5 < J < 2. At J = x loc 
obeys Curie's law. For J ^ but J < J c , we expect x loc 
to saturate at some cross-over temperature T* < J due 
to the coupling of the impurity spin to the planar mo- 
ments within the AFM LRO state. In agreement with 
this, Fig. signals a departure from x loc °t T~ x for 
T w 0.2 at J = 0.5, i.e. in the AFM LRO state. Sim- 
ilarly, for J > J c we expect a Pauli-like saturation of 
X loc for T <T* with T* < J due to the local dimer for- 
mation between the impurity spins and the planar sites. 
This can also be seen in Fig. for J = 2. The in- 
teresting point of Fig. however, is that it suggests a 
cross-over from the high-temperature Curie behavior to 
a region of power-law behavior x loc <x T~ a with an ex- 
ponent a different from unity in the vicinity of the QCP. 
Future QMC analysis should focus on additional data in 
the thermodynamic limit at T < 0.05 to elaborate on 
this observation. From Fig. we extract a « 0.20(5) 
at J = 1.40. The error on this exponent is rather large, 
due to the error in determining the QCP and due to the 



4 




Figure 6: Solid(dashed)[dashed-dotted]: ratio R of the total 
staggered structure factor and susceptibility vs. temperature 
for 0.5 < J < 2.0 and L = 24(50) [100]. In the classical 
(renormalized classical) regime, i.e. T — > oo (T — > 0, J < J c ), 
R is expected to be 1, In the gapped state Rl/T as T — * 0. In 
the quantum critical regime, roughly sketched by the shaded 
region, R differs from unity, however approaching a constant 
as T — > 0. Size of statistical errors is given by vertical bars. 

temperature range of only one decade to fit to. Regard- 
ing finite-site effects, the situation for x' oc is similar to 
that for Xu- As shown in the inset of Fig. EJ in the vicin- 
ity of the QCP the thermodynamic limit is reached for 
L > 24(50) if T > 0.1(0.05). Similar effects are expected 
at the lowest temperature T = 0.05 for J = 0.5 and 2.0 
and have not been considered. 

Finally we analyze the temperature dependence of the 
ratio of the total staggered structure factor to the to- 
tal staggered susceptibility of Eqn. Q. This is shown 
in Fig. El for 0.05 < T < 10 and 0.5 < J < 2. Non- 
linear error propagation of the QMC data through Eqn. 
© leads to substantially larger statistical errors on R as 
compared to the remaining quantities evaluated in this 
work. Several properties of R can be realized based on 
general grounds. First, for T ^ max {J, 1}, i.e. in the 
classical regime, \a (q) = Sa (q) /T for any wave vector 
q and therefore R — > 1. This behavior of R is obeyed for 
all values of J displayed in Fig. El Next, we note that 
the zero-temperature limit of the ratio Sa (Q) Ixa (Q) 



will be a T-independent constant whenever the system 
has no LRO at the wave vector Q and is gapped. This is 
true for any finite system, where Sa (Q) and xa (Q) will 
both saturate at finite values as T — > 0. It is also true in 
the thermodynamic limit where Sa (Q) and xa (Q) will 
both be exponentially activated to leading order. There- 
fore R oc l/T in the thermodynamic limit in the quantum 
disordered regime, i.e. for J > J c , which is consistent 
with the increase of R in Fig. El for J = 2. In addi- 
tion to this, R oc l/T for any other value of J below a 
characteristic temperature set by finite size gaps. This is 
particularly obvious for J w J c where strong finite size 
effects occur for T < 0.1. This effect also sets the mag- 
nitude of R for J » J c in Fig. Q] Finally, in the AFM 
LRO phase the systems allows for a classical description 
in terms of the order parameter modes leading to a renor- 
malized classical regime for which R (T — > 0) = 1 again 
US El- This is consistent with the behavior for J = 0.5 
in Fig. Eland with R in Fig. From a comparison of R 
at fixed T and identical J in panel a) and c) of the latter 
figure one can also deduce that the small difference be- 
tween R and unity in the renormalized classical regime 
decreases upon increase of L. Lowering the temperature 
from the classical to the renormalized classical regime, 
one crosses the quantum critical regime in which R > 1 
due to quantum fluctuations 0, |23| - For J < J c this 
regime has a finite extend in temperature only. Close 
to the QCP however, i.e. for J = 1.4 and 1.39, Fig. El 
strongly suggests that R approaches a temperature in- 
dependent constant R w 1.10(1) as T — > 0. Analysis of 
the NLcrM has resulted in R = 1.09 |S1 E^j. In turn 
the quantum critical regime starts at T ~ 1 and extends 
down to T — at the QCP. Unfortunately, in this regime, 
R is is very sensitive to the system size. This will be the 
issue of future QMC studies[18]. 
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